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1 Introduction 



Bethe ansatz method has proved to be the most powerful and (probably) unified method to 
construct the common eigenvectors of commuting families of operators (usually called trans- 
fer matrices) in two-dimensional integrable models pUEUHj. Faddeev et al pQ reformulated the 
Bethe ansatz method in a representation theory form: transfer matrices are associated with 
representations of certain algebras with quadratic relations (now called quantum groups). 
The eigenvectors are constructed by acting certain algebra elements on the " highest weight 
vectors" . This type of Bethe ansatz is known as the algebraic Bethe ansatz. 

Whereas this construction has been very successful in rational and trigonometric inte- 
grable models (HE], its extension to elliptic models had been problematic due to the fact that 
for the underlying algebras -Z n Sklyanin algebras (HI El the highest weight representations 
were not properly defined. Therefore, the algebraic Bethe ansatz method had not be applied 
to the elliptic integrable models directly. 

Recently, a definition of elliptic quantum groups E TjV (q) associated with any simple clas- 
sical Lie algebra g was given [S]. The highest weight representations of the elliptic quantum 
groups (cf . Z n Sklyanin algebras) are now well-defined [HI E] . This enabled Felder et al [TU] 
to apply successfully the algebraic Bethe ansatz method for constructing the eigenvectors of 
the transfer matrices of the elliptic integrable models associated with modules over E TjT] (sl2), 
and Billey ^T] to apply the algebraic nested Bethe ansatz method for constructing the eigen- 
values of the transfer matrices associated with a special module over E T ^(sl n ) (see section 
5). 

In this paper, we will extend the above construction of Bethe ansatz method further 
to the elliptic quantum group E T>rj (sl n ) on a generic higer spin module W = V^<; 1 )(z 1 ) £g> 
V A (i 2 )(z 2 ) <%> • • ■ <S> V A (i m )(z m ). After briefly reviewing the definition of the elliptic quantum 
group E TtV (sl n ), we study one parameter series of highest weight representations (higher spin 
representations) of E T}V (sl n ), and introduce the associated operator algebra and the transfer 
matrices corresponding to the E TjV (sl n ) -module in section 2. In section 3, we describe the 
algebraic Bethe ansatz for E T}V (sl n ). Finally, we give som e applications of our construction to 
the Z n Belavin models in section 4, and the elliptic Ruijsenaars-Schneider model associated 
to A n _i root system in section 5. 
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2 The elliptic quantum group and its modules associ- 
ated to A n _i 

2.1 The elliptic quantum group associated to A n -\ 

We first review the definition of the elliptic quantum group E T)V (sl n ) associated to A n _i 
[Hj. Let {ei | % = 1, 2, • • • , n} be the orthonormal basis of the vector space C n such that 
(ej, €j) = Sij. The simple roots are {on = — | z = 1, • • • ,n — 1} and the funda- 

mental weights {Aj | i — 1, • • • , n — 1} satisfying (Aj, a,-) = 5y are given by 

i . n 

jt=i fc=i 

Set 

i = ei — e, e= — z = l,---,n, then ^^ = 0. (2.1) 

n fc=l i=l 

For each dominant weight A = Y2a=i a *^-« > a * ^ there exists an irreducible highest 

weight finite-dimensional representation Va of with the highest vector |A). For example 

the fundamental vector representation is . In this paper, we only consider the symmetric 

i 

, * v 

tensor-product representation of Va 1 <8> Va 1 • • ■ (g> Vx ; (or, the higher spin-l representation of 
A ra _i), namely, the one parameter series of highest weight representations V^ii), with 

A« = /A 1; / G Z and I > 0. (2.2) 

This corresponds to the Young diagram I I I I • • • □. 

Let f) be the Cartan subalgebra of A n _i and f)* be its dual. A finite dimensional diagonal- 
izable f)-module is a complex finite dimensional vector space W with a weight decomposition 
W = ®tj,£f,*W\p], so that f) acts on W[/j] by xv = [i(x)v, [x G f), v G Wf//]). For ex- 
ample, the fundamental vector representation = C n , the non-zero weight spaces are 
W[i] =Cq, i = l,-..,n. 

Let us fix r such that Im(r) > and a generic complex number 77. For convenience, we 
introduce another parameter w = nn related to n. Let us introduce the following elliptic 
functions 



a[u) 



1 i 

2 
1 



(u,t), 6^\u) = 6 



1 _ 2 

2 n 
1 



Knr), (2.3) 
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a 
b 



(u, t) = ^ ex P { v^Tr [(m + afr + 2(m + a) (u + 6)] } . (2.4) 



m=— oo 



For a generic complex weight £ G J^™ =1 ^A*, we introduce a weight lattice with a shift 
by £, P^ as follows 



P^ = £ + niii, rrii G Z. 

For A G P 5 , define 

Aj=(A,ej), Ay = Aj - Xj = (A, e» - ej), i,j = !,-••, n. 



(2.5) 



(2.6) 



In this paper, we restrict A G Pg so that the inverse of the matrix S(z, A) in (J4.7)) does exist. 
Let R(z, A) G End(C n ® C n ) be the R-matrix given by 



A) = ^ A)£« ® + ^ A)^ ® E j3 + A)^ ® E i3 ) , 



i=l 



in which E i3 is the matrix with elements {E i3 ) l k = $jk$u- The coefficient functions are 



i$(*,A) = l, i$J(*,A) 



cr(z)cr(AjjW; + w) 
cr(,2 + w)a(Xijw) ' 



cr(w;)cr(,z + AjjU>) 



(2.7) 

(2.8) 
(2.9) 



a{z + w)a(Xijw) ' 

and Ay is defined in ()2.6|) . The R-matrix satisfies the dynamical (modified) quantum Yang- 
Baxter equation 



Ri2{zi - 22, A - h®)R 13 (zi - z 3 , X)R 23 {z 2 - 2 3 , A - h«) 
= #23(^2 - z 3 , X)R 13 (z 1 - z 3 , X - h ( - 2) )R 12 (z 1 - z 2 , A), 



with the initial condition 



i$(0,A) = 5l<5* 



(2.10) 



(2.11) 



We adopt the notation: Ru(z, A — h^) acts on a tensor v\ v 2 ® v 3 as A — fj,) ® Id if 

A representation of the elliptic quantum group E TjV (sl n ) (an £' T ^(s/ ri )-module) is by def- 
inition a pair (V4 7 , L) where J4 7 is a diagonalizable f)-module and L(z, A) is a meromorphic 
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function of A and the spectral parameter z G C, with values in End^C™ ® W) (the endo- 
morphism commuting with the action of fj). It obeys the so-called U RLL" relation 

R 12 ( Zl -z 2 ,X- h^)L 13 ( Zl ,X)L 23 (z 2 , A - fcW) 

= L 23 (z 2l X)L 13 (z u A - hW)R 12 ( Zl - z 2 , A), (2.12) 

where the first and second space are auxiliary spaces (C n ) and the third space plays the role 
of quantum space (W). The total weight conservation condition for the L-operator reads 

[h^ + h^, L 13 (z,X)} = 0. 

In terms of the elements of the L-operator defined by 

n 

L(z,X) (ei®u) = ^2e j ®li(z,X)v, veW, (2.13) 

3=1 

the above condition can be expressed equivalently as 

f(h)L((z, A) = L*(z, X)f(h + (2.14) 

in which f(h) is any meromorphic function of h and h measures the weight of the quantum 
space (W). 

2.2 Modules over E TyV (sl n ) and the associated operator algebra 

The basic example of an £' Ti?? (s/ n )-module is (C n ,L) with L(z,X) = R(z — zi,X), which 
is called the fundamental vector representation Va.i(^i) with the evaluation point Z\. It is 
obvious that ll RLL n relation is satisfied as a consequence of the dynamical Yang-Baxter 
equation ()2.10|) . Other modules can be obtained by taking tensor products: if (Wl ; L^) and 
(W 2 , L^) are E TiV (sl n ) -modules, where acts on (C n g> W s ), then also (Wi ® W 2 , L) with 

L(z,X) = L {1 \z,X-h {2) )L {2 \z,X) acting on C n ® W x ® W 2 . (2.15) 

An £' T ^(s/ n )-submodule of an i? r T? (s/ n )-module (W,L) is a pair (Wl,Li) where W\ is an 
f)-submodule of W such that C™ ® W\ is invariant under the action of all the L(z, A), and 
Lx(z, X) is the restriction to this invariant subspace. Namely, the Z? Ti ^(sZ n )-subrnodules are 
E TjV (s/ n )-modules. 
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Using the fusion rule of E TtV (sl n ) (|2.15|) one can construct the symmetric E TjV (sl n )- 
submodule of /-tensors of fundamental vector representations: 

VaV)( z i) = symmetric subspace of Vai(^i) <8> Vai( z i — w) ® • • ■ <S> V\i( z i — (I — ^) w ), 

where is defined by (|2.2jl . We call such an £ l Ti7) (sZ n )-module the higher spin-l represen- 
tation with evaluation point z\. These series of representations in the case of Z n Sklyanin 
algebra have been studied in jB] for n = 2 case and in ^21 E31 for generic n case. From direct 
calculation, we find that the E TV (sl n )-modu\e V^n(z) is an irreducible highest weight mod- 
ule of E T , n {sl n ) with the highest vector |A W ) G W[lAt] = C|A W ). It satisfies the following 
highest weight conditions 

l;(«,a)|a<')) = |aW), l<( u ,a)|a('))=o, » = 2,..., n , (2.16) 

L- (M ,A)|A(0> = ^ ^-^ A ^ + to ) |A(0>, iJ = 2,-,n, (2.17) 
J J a{u — z + iw) a{AiiW) 

/(/i)|AW> = /(/i)|AW>, (2.18) 

where /(/i) is any meromorphic function of h, which measures the weight of the quantum 
space W. 

For any i? Tfl (s/ n )-module, as in jTU] one can define an associated operator algebra of 
difference operators on the space Fun{W) of meromorphic functions of A G with values 
in W . The algebra is generated by h and the operators L(z) G End(C n <S> Fun(W)) acting 
as 

n 

L(z)(e t ® /)(A) = ® L fc, X )f( X ~ *)■ ( 2 - 19 ) 

One can derive the following exchange relation of the difference operator L(z) from the 
"RLL" relation (|2.12|) . the weight conservation condition of Ll(z,X) ()2.14|) and the fact 
[hW + h®, Ru(z, A)] = 0, 

#12(21 - 22, A - /i)Ia 3 (2i)X 23 (2:2) = L 23 (z 2 )L 13 (z 1 )R 12 {zi - z 2 , A), (2.20) 
f(h)L{(z) = Li(z)f(h + i-j), (2.21) 

where /(/i) is any meromorphic function of h and h measures the weight of the quantum 
space W. 
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The transfer matrices associated to an E T)V (sl n )-modu\e (W,L) [10J is a difference op- 
erator acting on the space _Fim(W)[0] of meromorphic functions of A with values in the 
zero-weight space of W. It is defined by 

n n 

T{u)f{\) = = I>i(«.A)/(A-<). ( 2 - 22 ) 

i=i i=i 

The exchange relation of //-operators (|2.20j) and (|2.21|) imply that, for any £^ T J? (s/ n )-module, 
the above transfer matrices preserve the space H = Fun(W)[0] of functions with values in 
the zero weight space VF[0]. Moreover, they commute pairwise on H: 



[T(u) , T(v)} 



= 0. 

H 



In this paper we will study the tensor product £y i7? (s/ n )-module W = V a q 1 )(zi) <g> 
V A (i 2 )(z 2 ) <8> ■ ■ ■ <8> V A (i m )(z m ). With the generic evaluation points {zi}, the module is an irre- 
ducible highest weight £ Tir ,(s/ n )-module [EJE]. Let A = H h A (u) , then W[A] = C|A) 

with 

|A) = \A {h \ A {h \ • • • , A (u) ) = |A (Zl) ) ® |A fe) ) ® • • • ® |A (im) ). (2.23) 

The vector |A), viewed as a constant function in Fim(W^), obeys the following highest weight 
conditions: 

L\(z)\A) = A(z,X)\A), L\(z)\A) = 0, i = 2,---,n, 

L){z)\A) = 5)D l {z,\)\A), i,j = 2,---,n, f(h)\A) = /(M)|A). 

The highest weight functions read 

A( \\ i nf \\ )YT^(z-pk)\(r(X ll w + Nw) . 

A(z,X) = l, D i (z,X) = <\\— -} r , i = 2,---,n, (2.24) 

where 

m 

Pk = z k , q k = z k -l k w, N = ^l k , fc = l,---,m. (2.25) 

fe=i 

3 Algebraic Bet he ansatz for E TjV (sl n ) 

In this section we fix a highest weight £' TiI/ (s/ n )-module W of weight A, the functions 
A(z, X),Di(z, A) dOD) , with the highest vector |A). We assume that AT = ^™ =1 l k = n x I 



with I being an integer, so that the zero- weight space W[0] can be non-trivial and so that 
the algebraic Bethe ansatz method can be applied as in [T4*l \W\ \W\ HUj. 
Let us adopt the standard notation for convenience: 

A(u)=L{(u), B i (u) = Ll(u), i = 2,...,ra, (3.1) 
b(u)=L[(u), Vi(u)=Li{u), i,j = 2,---,n. (3.2) 

The transfer matrices T(u) become 

n 

T(u) = A(u) + J2vi(u). (3.3) 

i=2 

Any non-zero vector \Q) G -F-un(W)[A] is of form |0) = g(X)\A), for some meromorphic 
function g ^ 0. When N = n x I, the weight A can be written in the form 

n-1 

A = n/Ai = iJ3(ei-e Jfe+ i). (3.4) 

k=l 

Noting (|2.21j) . the zero- weight vector space is spanned by the vectors of the following form 

Bw 1 (vjVi) B iw 1 -i( u JVi-i) ' • (3.5) 

where Nx = (n — 1) x I and among the indices {i& | k — 1, • - • , iVi}, the number of = j, 
denoted by should be 

#(;) = /, j = 2, ■ ■ ■ ,n. (3.6) 
The above states ()3.5j) actually belong to the zero-weight space W[0], because 

iVi 

= ^kjs^^k^o • • ■B il {v 1 )f(h + J2ik - Nxi)\n) 

jt=i 

n 

= /(^« - (n - 1)H + nri^KJS^^-i) • ■■B il (v 1 )\ty 

i=2 

for any meromorphic function /. We will seek the common eigenvectors of the transfer 
matrices T(u) in the form 

|A; K» = f*iA.~M(a ; KDB^KJB^^K^x) • • ■ B^n), (3.7) 
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with the restriction condition ()3.6|) and the parameters {vk} will be specified later by the 
Bethe ansatz equations (|3.46j) . We adopt here the convention that the repeated indices imply 
summation over 2, 3, • • • , n, and the notation that 

v k = v ( °\ k= 1,2,- (3.8) 

For convenience, let us introduce the following set of integers 

Ni = (n-i)xl, z = l,2,---,n-l, (3.9) 

and (n — 1) complex parameters {a.® |z = 1, • ■ ■ , n — 1} and set = — ^2=1 a ^ ^° specify 
quasi-vacuums of each step of the nested Bethe ansatz (see below). Associated with {a^}, 
set 

a« = - Ja' i+1 '- Lt=i+1 L i = 0,.-.,n-2. (3.10) 

(n — z — 1) I n — i 

From the exchange relations (|2.20|) and (|2.21|) . one can derive the commutation relations 
among A(u), T>f(u) and Bi(u) = 2, • • • ,n) (for details, see Appendix A). The relevant 
commutation relations are 

A(u)B t (v) = r(v - u, X il )B i (v)A(u) + s(u - v, X u )Bi(u)A(v), (3.11) 

n 

Vi(u)B l (v) = r(u-v, \ n -h n ) Rf l (u-v,X)B^v)Vi(u) 

a,(3=2 

- s{u - v, Ay - h^Bii^Vjiv), (3.12) 

n 

Bi{u)Bj(v) = R fM ~ v > X)Bf,{v)B a (u), (3.13) 

a,P=2 

f(h)A(u) = A(u)f(h), f(h)Bi(u) = Bi{u)f(h + 1-1), (3.14) 
f(h)Vi(u) = Vi(u)f(h + i-j), (3.15) 

where the function r, s are defined by 

, . a(u + w)a(uw) , . a(w)a(u + ww) _ 

r(u,u)= ' y , 8 (u,u) = -)-^ L ueC. 

a{u)a{pw + w) G\u)<j{yw — w) 

Because of the simple poles in the functions r, satnG Z + tZ, let us assume for convenience 
that all {v^ are distinct modulo Z + rZ and consider u at a generic point. Take 



|n) = 0(A)|A), 

9 



(3.16) 



as the so-called quasi- vacuum satisfying the following conditions 

A(u)\n) = 9{x ~ 1) \n), d{u)\n) = o, % = 2, • (3.17) 

9\ A ) 



<?(A) 

/(A)|J2) = /(JVi)|n>, Bi(«)|fi)^0, i = 2,...,n, (3.19) 



^( l i)|fi) = 5;AKA)^|fi), z,j = 2,---,n, (3.18 



where A), and are given in (|2.24jl and (|2.25jl . We will specify the function g(X) 

later (l3~23l 

Now, let us evaluate the action of A(u) on |A; {ffc})- Many terms will appear when we 
move A(u) from the left to the right of Bi(vk)'s. They can be classified into two types: wanted 
and unwanted terms. The wanted terms in A(u)\X; {vj~}) can be obtained by retaining the 
first term in the commutation relation (J3.11|) . The unwanted terms arising from the second 
term of (|3.11|) . have some B(vk) replaced by B(u). One unwanted term where B(vni) is 
replaced by B(u) can be obtained by using firstly the second terms of (J3.11j) . then repeatedly 
using the first term of (j3.ll)) . Thanks to the commutation relation (|3.13jl . one can easily 
obtain the other unwanted terms where B(vk) is replaced by B{u). So, we can find the action 
of A(u) on I A; {ffe}) 

A(u) F^-^{\- {v k })B iNi (vnJB^ (wjvx-i) ■ ■■B il (v 1 )\Q) 

= F^-^ (A - 1; {v k }) \ II - (A - E L Ci h ~ <*» 

x B iNi (v Nl ) ■ ■ ■ BM^^-in) 



'JVi-l Ni 



x B lNi (u)B tNi ^(v Nl ^) ■ ■ ■ Bil ( Vl )^p-\n) 
+ o.u.t., (3.20) 

where o.u.t. stands for the other unwanted terms. Due to the restriction condition ()3.6|) and 
Ni = (n — 1) x I, we have 

n , 1 j -pi o-(t> fe -m + w) [ J A or(Ajitu - Ztw + w) 



. fc=l s=fc+l ,/ k&=l J U=2 



L cr^-ii) J 1 A = A a(X jl w + w) 
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(3.21) 

Ni n 

Y^i k = = -li. (3.22) 

k=l i=2 

In order that the Bethe ansatz equation ()3.46j) and the eigenvalues (|3.31|) should be inde- 
pendent of A, one needs to choose (cf. [TP] ) 

.(A)=^-»n(n "X) h,) l- < 3 ' 23 > 

3=2 U=l K J ) 

The action of A(u) on |A; {i>fc}) becomes 

A(u) F h ^-^(X; {vkVBi^VNjBi^vn-i) ■ ■■B il (v 1 )\ty 

= ^ aWw ft " ( ^-^r ) Fn ' t2 ''''' tiVl(A - i; H^KJ • ••S il (t;i)|n) 

/ x , o-(A ijv + w) (i) ^'J-r 1 a(u fe - Ujvi + iw) 
+ s (u-v Nl ,Xu N ) tt » II — 

x /•'''•'•••• -' N : (A - 1; {M^HB^^-i) ■ 
+ o.ui.. (3.24) 

Next, we evaluate the action of T>\{u) on |A; {vk})'- 

Ni 

Y[ r{u - Vh ,(X-h- (JV X - Q - ei)) 
,fc=i 

m 



X 



-p-r <r(tt -pfc) I q~(Ajiw + (N - l)w) g(X-i + 11) 
L\*(u-q k )j a(XuW-lw) g(X + ll) 



k 



i j' ... ;,i 



x A) ^'-^(A-ijK}) 

l. ) 11,12, •••,ijv 1 

+ u.t., (3.25) 
where it.i. stands for the all unwanted terms. We have introduced the convenient notation 

il,i2 ' ' JJVl k=2 k=3 

x-.-Rj^(u-v Nl ;X), (3.26) 
11 



and all the indices take the values 2, • • • ,n. Noting the explicit form of the function g(X) 
given in (j3.23|) and the zero-weight of the state |A; {v k }), the action of V\{u) becomes 

Vi(u)\\;{v k }) 



e ™ 



( 21 W J rr <j(u -v k + w) \ J-rj a(u-p k ) 
11 a (u-v k ) j\lla(u-q k ) 



{' j' ... i' 



k\L^I(u,\)\ 1 F**»™*«i{\-i;{v k }) 

x B^vn^B^ ^vn^) ■■■B i , 1 {v 1 )\n) 

-en s(u-v Nl ,X u ) [[ — {{ 



+ o.u.t.. (3.27) 

Keeping the initial condition of R-matrix (j2.11|) and the notation ()3.26|) in mind, we can find 
the action of the transfer matrices on the state |A; {v k }): 

T(u)\X; K» = [A(u) + J>K«) ) |A; {v k }) 



i=2 



fc= 
iVi 



+ e 41 w I T-T g(M ~ ffc + W) I I -pr 0-(M-p fc ) 

x 8^(^1)^-1 ( v *i-0 • ■■B l[ (v 1 )\n) 

i \ , (T (^jv 1 iW + w ) k» a (i) m ^t-t 1 cr(v k - v Nl + w) 
+ s(u-v Nl ,X u ) — e » M — 

^i Nl i w ) fe = i <r(v k -v Nl ) 

x F u,i2,- -,% l(A _ i- { Ufc })s ijVi ( M )5 ijVi _ 1 ( UjVl _ 1 ) • • -B^vJin) 

E<*^ w , . o-(Aj!W + w) J - v fc + to) 1 J -A- a(v Nl - p k ) 

% ._r s{u - vN - Xu) ^x-^r{ji a( VNl - Vk) }{n 

x^I^kj};;;::;:;:;; ^••■^(a ; w}) 
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+ o.u.t., (3.28) 
where we have introduced the reduced transfer matrices T^(u) 

n 

(3.29) 



i=2 



Thanks to the commutation relations ()3.13|) . one can obtain the explicit expressions of all the 
o.u.t. The equation 1)3.28)1 tells that the state |A; {vk}) is not an eigenvector of the transfer 
matrices T(u) unless F's are the eigenvectors of the reduced transfer matrices T^\u). The 
condition that the third and fourth terms in the above equation should cancel each other 
and also the all o.u.t. terms vanish, will give a restriction on the N\ parameters {vk}, the 
so-called Bethe ansatz equations. Hence we arrive at the final results: 

T(u)\X; K» = t(u; K})|A; {v k }). (3.30) 

The eigenvalue reads 

i(«;W) = e "-)°" (ft Jfa ~ U + W) 



N 



L a(v k - u) 



(3.31) 



in which 



rW( w )F*i^'-^i(A; H» = e-^r"t (1) ( M ; (A; {v fe }), (3.32) 

^^.--^(a _ i ; { Vk }) = ^F ii '^'- > ^i(A; (3.33) 



and {^fc | fc = 1, 2, • • • , iVi} satisfy 



g—naw 



^-, + ») J n ^U (3 . 34) 



The parameters k — 1, 2, • • • , iV 2 } will be specified later by the Bethe ansatz equations 

dug} . 

The diagonalization of the transfer matrices T(u) is now reduced to diagonalization of the 
reduced transfer matrices T^(u) in ()3.32)1 . The explicit expression of T^'(u) given in ()3.29)) 
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implies that T^'(u) can be considered as the transfer matrices of an £' T ^(s/ n „i)-module 
(or reduced space): Ni tensor product of fundamental representations of E TtT] (sl n _i) with 
evaluation points {vk}- We can use the same method to find the eigenvalue of T^(u) as we 
have done for the diagonalization of T(u). Similarly to (j3.1|) — ()3.2}) . we introduce 

A^\u) = {L^}l(u), B?\u) = {L^{u), < = 3,...,n > (3.35) 
Ci 1 \u) = {L^y 2 (u), vW(u) = {L^y i (u) 1 i,j = 3,...,n. (3.36) 

Then the reduced transfer matrices ()3.29|) can be rewritten 

n 

TW(u) = A (1 \u) + ^V il ^{u). (3.37) 

i=3 

We seek the common eigenvectors of the reduced transfer matrices T^'(u) in an analogous 
form to ()3.7j) : 

|A;{4 1} }) (1) = £ i^ 2 '"^ 

(3.38) 

with the restriction of zero- weight conditions similar to ()3.6|) . The quasi- vacuum is 
the corresponding highest weight vector in the reduced space We can also derive the 

relevant commutation relations among A^(u), T>^\ and B^(u) = 3,4, • • ■ ,n) similar 
to (|3.11|) — (|3.15|) . Using the same method as we have done for the diagonalization of T(u), 
we obtain the eigenvalue of the reduced transfer matrices T^(u): 

T {1) (u)\X; = {^ 1} })|A; (3.39) 

Though very complicated, the coefficients F %1 ' %2 ' '"> %N i (A; {v^}) could be extracted from the 
eigenvectors |A; {v\, })^, in principle. The eigenvalue ^^(w; {f^ }) is given by 

f ^2 / (1) | \ 



(3.40) 
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(2) 

The parameters {v k \ k = 1, 2, • • ■ , A^} will be specified later by the Bethe ansatz equations 
(|3.46|) . The function t^ 2 \u; {v^}) is the eigenvalue of the second reduced transfer matrices 
T^{u) 

T (2) (M)j p(D^-^ 2(A . {v (i) }) = e -BfA^(u; {v^})F^'-'^(X; {4 X) }),(3.41) 

and T( 2 >(u) is given by ()3.26|) and ()3.29|) with all indices taking values over 3,---,n and 
depending on {v k } instead of {v k }. Hence, the diagonalization of T(u) is further reduced 
to the diagonalization of T^ 2 \u) in (j3.41j) . This is so-called nested Bethe ansatz. Repeating 
the above procedure further n — 2 times, one can reduce to the last reduced transfer matrices 
T^ n ~ l ^{u) which is trivial to get the eigenvalues. At the same time we need to introduce the 
n ^ n ~ 1 - > / parameters {{v k \ k = 1, 2, ■ • • , N i+ i}, i — 0, 1, • • - , n — 2} to specify the eigenvectors 
of the corresponding reduced transfer matrices T^\u) (including the original one T(u) = 
T^(u)), and n arbitrary complex numbers {a^i = 1- • - ,??,} to specify the quasi-vacuum 
of each step as in 1)3.16)1 and 1)3.23)1 . Finally, we obtain all the eigenvalues of the 
reduced transfer matrices T^(u) with the eigenvalue t^(u; {v k }) in a recurrence form 

t il) (u;{v^}) = e^ l+1 - n ^ )w 

{Ni+i , (i) \} ( Ni , (i) s ^ 

g( + > | fjjf^Jl ,<. +1)(M; {„<«.}), 
fc= l <J\u-v\>) J \ k=1 o(u-q\ r )) 

i = 0,l,---,n-2, (3.42) 
t^\u- {^ 0) }) = t{u- {v k }), &- l \u) = 1, (3.43) 

where a«, i — 0, 1, • • • , n — 2 are given by ()3.1(Jj) . a = a(°), N = m and 

Pk )= Pk = z kl qf> = q k = z k - l k w, k = 1,2, - ■ ■ ,m, (3.44) 
vf = v { k ' 1 \ qf=vt l) -w, i = l,2,...,n-2, fc = 1, 2, • • • , N t . (3.45) 

The {{i>jt }} satisfy the following Bethe ansatz equations 

We conclude this section with some remarks on analytic properties of the functions 
t®(u; {v k })- By construction, the eigenvalue functions t^(u; {v k }) of the transfer ma- 
trices should not be singular at m = uj^ (modulo Z + rZ) with < i < n — 2, l</c< A^+i. 
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On the other hand, the Bethe ansatz equations (|3.34|) and ()3.46|) are derived from the require- 
ment that the unwanted terms should vanish. It is interesting to note that these constraints 
could be understood from a different point of view: from equation (|3.31jh we know that 
u = = Vk (modulo Z + rZ) is a possible simple pole position of t(u; {vk})- However, 
the constraints on {v^}, the Bethe ansatz equations (|3.34|) simply tell that the residue of 
t(u; {vk}) at u = Vk (modulo Z + rZ) is vanishing. Hence, t(u; {vk}) is analytic at u — Vk 
(modulo Z + rZ). Similarly, the Bethe ansatz equations ()3.46|) ensure the analyticity of 
t^(u]{vf}) at all u = Vj? (modulo Z + rZ). Therefore, the eigenvalues of transfer matrices 
T(u) are analytic functions of u. 

4 The r L n Belavin model 

We shall show in this section how to obtain the eigenvectors and the corresponding eigen- 
values of the transfer matrices for Z ra Belavin model with periodic boundary condition from 
our results in section 3. 

The Z n Belavin model J7j is based on the following R-matrix 

Rb(u)= J2 R %( u ) E i*® E ir 

i,j,k,l 



The coefficient functions are [T5] 

h(u)cr(w)e^-^ (u+w) 



Rf.(u) = <J v(u+ w )9(>-V(w)8(k-j)(u) rfi+j-k + l mod n 
otherwise 



where 



u 



h( u ) = , ; \ (4.3) 

The R-matrix satisfies the quantum Yang-Baxter equation 

R 12 (ui - u 2 )R 13 (u 1 - u 3 )R 23 (u 2 - u 3 ) = R 23 (u 2 - u 3 )R 13 (ui - u 3 )R 12 (ui - u 2 ), (4.4) 

and the initial condition: Rij(0) = SjSj. For generic evaluation points z%, ■ ■ ■ , z m , one then 
defines the commuting transfer matrices 

T B (u) = tr L B {u), (4.5) 

L B (u) = {R B }oi(u - Zx) ■ ■ ■ {RB}om(u - Zm), (4.6) 
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where the L-operator Lb(u) acts on C n <g> (C n ) m and the transfer matrices Tb(u) act on 

Using the intertwiner introduced by Jimbo et al [T%] . we can define the matrix S(z, A) 
with the elements £^(2, A) = ^ ^^(X^) • When A G Pg, the inverse matrix of S(z, A) exists. 
Then, one finds the following twisting relation holds: 

{R B }i2(u 1 - u 2 )S( Ul , \) {1) S(u 2 , A - hU)W 

= S(u 2 , A) (2) S(nx, A - h^)^R 12 ( Ul - u 2 , A), (4.7) 

in which R(u,X) is defined by ()2.7|) . ()2.8|) and ()2.9|) . The modified quantum Yang-Baxter 
equation (j2.10J) and the quantum Yang-Baxter equation (j4.4j) are equivalent to each other 
due to the above twisting relation ()4.7|) . Let us introduce 

m 

S m ( Zl , ■■■ ,z m ;X) — S(z m , X)^S(z m ^, A - h^)^ ■ ■ ■ S(z 1} A - ^ h^f\ (4.8) 

k=1 

For the special case of li = l 2 = ■ ■ ■ = l m = 1, we find that the L-operator L(u,X) 
corresponding to the £' T)r/ (s/ n )-module W is equivalent to the L-operator Lb{u) of Z n Belavin 
model ()4.6|) by the twisting relation 

S m (zi, ■■ ■ ,z m ; X)S(u, X - /i) (0) L(w, A) 

= L B {u)S(u, X)^S m ( Zl , •■■ ,z m ;X- h®). (4.9) 

In this formula S m (zi, • • ■ , z m ; A) acts on the factors from 1 to m, and h = + ■ • • + h^ m \ 
Acting the both sides of (|4.9|) on a state ej <8> v with v G W[0], replacing A by A + i, and 
noting (J2.14|) . we find 

n 

^2 (s(u,X + j)e^j ® (s m (z u - ■■ ,z m \ X + i)Ll(u, X + i)v) 

= L B {u) | (S(u, X + i)eij <g> (S m {zi, ■■■,z m ; A)u)} . (4.10) 
Therefore, one can derive 

n 

T B {u)S m ( ; A)|h — »S m ( X + i)L}(u,X + i)\ H , (4.11) 

i=l 

where H = Fun(W)[0] is preserved by the transfer matrices. 
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Let / i — > J f(X) be a linear functional on the space T of functions of A G C™, such that 
J /(A + i) = J /(A), % = 1, • • ■ , n, for all / G T. Extend J to vector-valued functional by 
acting component-wise. Then for each eigenfunction \1/(A) of T(u) with eigenvalue t(u), the 
transfer matrices of Z n Belavin model Tb{u) act on the vector J S m (zi, ■ ■ ■ , z m 

T B (u) I S m (zi, ■•■,z n ; A)tf (A) = ^ / S^zi, • • • , z m ; A + i)Ll(u, A + (A) 
J i=i ^ 

n » 



S m {z 1 ,---,z m ;X)T{u)^{X)=t{u) J S m (z x , ■ ■ ■ , z m ; A)tf(A). (4.12) 

Namely, the vector J S m (zi, • • • , z m ; A)\1/(A) is an eigenvector of T B {u) with the same eigen- 
value. Fortunately, such a linear functional on T has been given in QU] for XYZ spin 
chain (Z^-Belavin model) and its higher spin generalizations [TSj, for the generic Z n Belavin 
model ^3]. And our results recover those of JUj for n = 2 case and coincide with those of 
[To] for Z n Belavin model. 

As in [T^], one can construct the higher spin generalization of Z n Belavin model by 
putting at each site a local L-operator 



L k {u- l k ) = V a ° [U ** +J I a <g> a {5 Q } , / fc G Z+, <S Q G (4.13) 
z — ' a, * — > 



• n ■ 



where pi k is the spin-4 representation of Z„ Sklyanin algebra SK with dimension ^zj^r 
[T2*] . The corresponding transfer matrices of the model are given by 

T B (u; {l k }) = tr {L±{u- h)L 2 (u; l 2 ) ■ ■ ■ L m (u; l m )} . (4.14) 

Then, our construction in Sect. 3 actually results in the eigenvalues formulation ()3.31|) of the 
transfer matrices (|4.14|) . with the Bethe ansatz equations ()3.46|) . In particular, our result 
recovers that of 16 j when n — 2. 

5 The elliptic A n -\ Ruijsenaars- Schneider model 

If we take £' r ^(si n )-module W for the special case 1 

W = V A(h) (z 1 ), A {h) = ZiAi = A^Ai = n x IA X , (5.1) 

1 Billey used algebraic Bethe ansatz method over the -E Tjr) (si„)-module W — V A (i)(z 1 ) ® • • - V A (i)(zi — 
(nl — l)w) to obtain eigenvalues of elliptic N-body Ruijsenaars operator mj 
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the zero-weight space of this module is one-dimensional and the associated transfer matrices 
can be given by 

= g (« - zi + Iw) M 
a(u — zi + nlw) 
The operator M is independent of u and is given by 

Here {1^} are elementary difference operators: Tif(X) = f(X — i). In fact, for a special 
choice of the parameters [T^], this difference operator M is the Hamiltonian of elliptic A n ^\ 
type Ruijsenaars-Schneider model [20] with the special coupling constant 7 = Iw, up to 
conjugation by a function |3TJ. Now, we consider the spectrum of M. The results of 
Sect. 3, enable us to obtain the spectrum of the Hamiltonian of the elliptic Ruijsenaars- 
Schneider model as well as the eigenf unctions, in terms of the associated transfer matrices 
()5.2j) . Since M is independent of u, we can evaluate the eigenvalue of T{u) at a special 
value of u, u = z\. Then the expression of the eigenvalue t{u\ {vk}) simplifies drastically, 
for the second term in the right hand side of ()3.31|) (the one depending on the eigenvalue of 
the reduced transfer matrices t^(u; {v k ^})) vanishes because u — pf' = in the case ()5.8|) . 
Finally, we obtain the eigenvalues of M: 



a(lw) ) a{v k -zi) 



n v r v ■ < 5 ' 4 > 

J- J- rr in — 7, 



where {^} and {{fj^}} are determined by the nested Bethe ansatz equations 

I ^+1 / (i) 00 , \ "I f ^ / (*) I 

'-**"•" n g V" ) - n frf {<r>». 

-w) J U=i°"(^ 

(5.5) 

The functions tM (^; {-u^ }) appearing in (|5.5|) are defined by the following recurrence relations 

W+i / (i) 




1 00 \ 

fc =i -m) 

+ea ,.,„ f jj g (u-« t + fjj 1 i(>+1)(M; { „, +1)}) , 

i = l,---,n-2, (5.6) 
f(™- 1 )( M ) = l, (5.7) 
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where i = 0, 1, • • • , n — 2 are given by (|3.10|) . a = a^°\ N = 1 and 

pf^ = zi, = Z\ — nlw, (5.8) 
vf=vt 1 \ g« = 4 4 - 1} -«;, i = l,2,...,n-2, fc = 1, 2, • • • , iV,. (5.9) 

6 Conclusions 

We have studied some modules over the elliptic quantum group associated with A n _i type 
root systems. There are elliptic deformation of symmetric tensor product of the fundamen- 
tal vector representation of A n _i (we have called it higher spin representation). Using the 
nested Bethe ansatz method, we exactly diagonalize the commuting transfer matrices asso- 
ciated with the tensor product of the £ , T)7? (s/ ri )-modules with generic evaluation points. The 
eigenvalues of the transfer matrices and associated Bethe ansatz equation are given by (|3.31|) 
and ([3.46)1 . For the special case of n — 2, our result recovers that of Felder et al [TU] . 

We also take the applications of our construction to some integrable models. For li = 
l 2 — ■ ■ ■ — 1 case, our result gives the spectrum problem of Z n Belavin model with periodic 
boundary condition, which coincides with the result [THj. For the special case W = V A ( n i) (zi), 
or equivalently 1% — • ■ • — l n i — 1 and Zi = Zi — (i—l)w, our result gives the spectrum problem 
of the elliptic A n _i type Ruijsenaars operators, which is "relativistic" generalization [22] of 
the Calogero-Moser type integrable differential operators, with the special coupling constant 

7 = Iw, which coincides with the result jTIJ. Moreover, for the generic case, our result gives 
eigenvalues of the transfer matrices given by (j4.14J) associated with higher spin generalization 
of Z n Belavin model with periodic boundary condition, and it further enables us to construct 
the eigenvalues of all types of Ruijsenaars operators associated with the A n _i root system 
in the Bethe ansatz form [T9*] . 

Acknowledgements 

We thank professors A. Belavin, S. Odake and K.J. Shi for their useful discussion. This work 
is supported in part by Grant-in- Aid for Scientific Research from the Ministry of Education, 
Culture, Sports, Science and Technology, No. 12640261. W.-L. Yang is supported by the 
Japan Society for the Promotion of Science. 



20 



Appendix A: The commutation relations 

The starting point for deriving the commutation relations among A(u), and Bi(u) 

= 2, • • • , n) is the exchange relations ()2.20|) and (j2.21|) . We rewrite ()2.20|) in the com- 
ponent form 

n n 

-v,\- h)Li(u)L^(v) = L b f (v)L$(u)rt j, d (u - v, A). (A.l) 
i',j'=i i'j'=i 

For a = b = d = 1, c = i ^ 1, we obtain 

it-{(M — f , A — 1 — i) R\i\u — v, A — 1 — i) 

' A(u)A(v) - - V ' ^- Bi(v)A(u), (A.2) 



B\\(u-v,X) w w i$(u-v,A) 

by noting the definitions ()2.8|) and ([2.9)1 . The commutation relation (|3.11j) is a simple 
consequence of the above equation. 

Similarly, for a = j ^ 1, b — 1, c = « ^ 1, d = I ^ 1, we obtain 

t ( i- ?(M ~ - " r wm. (a.3) 

I ^l( M ~V,X-h) J ^(U -V,X-h) 

from which follows the commutation relation ()3.12j) . For a = b = 1, c = i ^ 1, d = j ^ 1, 
we obtain 

n n 

BiMBjiv) = Bp{v)B a {u)R«?{u -v,X) = J] R°f(u - v, X - a - $)B p (v)B a {u) 

a,f3=2 a,0=2 
n 

= R ff( u - v -> ^)Bp{v)B a {u), (A.4) 

a,/3=2 

which leads to the commutation relation (|3.13j) . The other commutation relations ()3.14j) - 
([3.15)1 are derived from (|2.21j) . 
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